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Effect of surface nanostructure on temperature programmed reaction spectroscopy: 
First-principles kinetic Monte Carlo simulations of CO oxidation at RuO2(110) 
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Using the catalytic CO oxidation at RuO2(110) as a showcase, we employ first-principles kinetic 
Monte Carlo simulations to illustrate the intricate efi'ects on temperature programmed reaction 
spectroscopy data brought about by the mere correlations between the locations of the active sites 
at a nanostructured surface. Even in the absence of lateral interactions, this nanostructure alone 
can cause inhomogeneities that cannot be grasped by prevalent mean-field data analysis procedures, 
which thus lead to wrong conclusions on the reactivity of the different surface species. 



PACS numbers: 



;.43.Bc, 68.43.De, 68.43.Vx, 82.65.-l-r 



Temperature programmed desorption (TPD) and re- 
action (TPR) spectroscopy is a frequently employed tool 
that provides unparalleled insight into the binding en- 
ergetics of adsorbates or reactants at solid surfaces by 
recording the amount of desorbing species while ramp- 
ing the substrate temperature [1[. The obtained experi- 
mental data does contain information about the spatial 
arrangement of the surface species before and during the 
ramp, but only in an indirect way. Prevalent experimen- 
tal data analysis procedures neglect such dependencies 
and rely almost entirely on a mean-field approach, arriv- 
ing at expressions considering only the averaged surface 
coverages 0i of species i [l[ . Obvious breakdowns of the 
mean-field assumption come from strong lateral interac- 
tions among the adsorbates or extrinsic defect sites. The 
intricate consequences on TPD and TPR spectra are well 
documented in the literature [2] and have made surface 
scientists broadly aware of the limitations of the preva- 
lent mean-field data analysis procedures. With the ex- 
plicit lattice structure not taken into account, another 
fundamental source of deficiency of mean-field theories 
are the neglected correlations between the locations of 
the active surface sites themselves, i.e. that the mere 
way how the active sites are arranged on the surface is 
not accounted for. Within the general move towards the 
atomic-scale characterization of ever more complex and 
nanostructured surfaces exhibiting a wide variety of dif- 
ferent surface sites this will become an increasingly im- 
portant factor also for the analysis of TPD and TPR 
spectra. We illustrate this here for the catalytic CO ox- 
idation at RuO2(110), which is a showcase for a system 
with only very small lateral interactions. Without appre- 
ciable thermodynamic driving force for segregation (and 
thus for inhomogeneities in the adlayer) one would gen- 
erally expect a mean-field picture to be valid. Instead 
we demonstrate that the neglect of the essentially one- 
dimensional trench like arrangement of the most active 
surface sites in a mean-field analysis of a set of experi- 
mental TPR data |3J for this system leads to qualitatively 
wrong conclusions about the reactivity of the various sur- 
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FIG. 1: Top view of the RuO2(110) surface showing the 
two prominent active sites (bridge and cus). Ru = light, 
large spheres, O — dark small spheres. When focusing on 
these two site types, the surface can be coarse-grained to the 
lattice model on the right, composed of alternating rows of 
bridge and cus sites. Atoms lying in deeper layers have been 
whitened in the top view for clarity. 



face species. 

Over the past years, RuO2(110) and the CO oxida- 
tion at this model catalyst surface have been extensively 
studied both experimentally, as well as by first-principles 
theory [4| . A wealth of data from ultrahigh vacuum ex- 
periments and density-functional theory (DFT) calcula- 
tions is complemented by detailed in situ measurements 
and first-principles kinetic Monte Carlo (kMC) simula- 
tions of the steady-state catalytic activity. The picture 
that is unanimously obtained is that the surface kinetics 
is predominantly taking place at two prominent active 
sites, the so-called coordinatively unsaturated (cus) and 
the bridge (br) site. At the rutile (110) surface these 
two sites are arranged in alternating rows as illustrated 
in Fig. [1] Of the corresponding two potentially active 
oxygen species, 0'^"'^ and O''', it is particularly the much 
weaker bound (by r^ 1.5 eV) cus oxygen that would be 
expected to be mostly responsible for the high catalytic 
activity of this model catalyst at high pressures. This 
general expectation derives from the established impor- 
tance of the oxygen- metal bond breaking step in catalytic 
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FIG. 2: (Color online) Total CO2 yield from surfaces initially 
prepared with varying 0"^"° coverage 9. In all surfaces all 
bridge sites are initially covered with O*"^ species, and the 
remaining (1 — ^) cus sites not covered by O'^"'' are occupied 
by CC^^^. The CO2 yield is given relative to the one obtained 
for the surface with zero 0"^"° coverage. Shown are the total 
simulated CO2 yield, as well as the contributions from the 
two dominant reaction mechanisms under these conditions: 
Qcus_^(-,Qcus g^jjj O'^'+CO'^"". The experimental data is taken 
from ref. [3|. 



cycles and the notion behind the weU-known Br0nsted- 
Evans-Polanyi type relationships {5]. In recent detailed 
first-principles statistical mechanics modeling this gen- 
eral expectation has been fully confirmed 6]: The cal- 
culated O'^^'^-fCO'^"'' barrier is with 0.9 eV significantly 
lower than the one of the competing O'''' -I- CO'^"*' 
tion (1.2 eV), and the prior process involving the O 
species was indeed found to completely dominate the 
simulated steady-state activity of this surface at tech- 
nologically relevant gas phase conditions. A recent set 
of TPR experiments came therefore as a particular sur- 
prise, since the observed dependence of the CO2 yield on 
surface coverage was put forward as evidence for a com- 
parable activity of the strongly-bound O*^' surface oxy- 
gen species |3| , thereby implying either a violation of the 
general Br0nsted-Evans-Polanyi relationships or errors in 
the relative first-principles barriers of the competing reac- 
tions that exceed 0.3 eV (and thereby amount to several 
orders of magnitude in the temperature-dependent rate 
constants) . 

Specifically, TPR spectra were recorded from surfaces 
in which initially all bridge sites were always fully covered 
with O*^', while the coverage of O'^"'' varied in the range 
Q < 9 < 0.8 monolayer (ML), where 1 ML corresponds 
to an occupation of all cus sites. The remaining free cus 
sites were then each time saturated with CO, so that the 
initially prepared surfaces contained an amount of IML 
0"°', e ML 0^="" and [1-9) ML C0'="^ Within the under- 
standing of a much more reactive 0'^"'^ species and thus a 



dominance of the O'^^'^-I-CO'^^'' reaction, the anticipated 
coverage dependence of the total CO2 yield in the TPR 
experiments in a mean-field picture is simply a parabolic 
9{1 — 9). What was measured instead is an almost linear 
(1 — 6') decrease of the total CO2 yield with increasing 
Qcus coverage as shown in Fig. [2] [3|. Such a behavior is 
much easier rationalized assuming the O*"^ species to be 
more reactive (and correspondingly the O'^'^-I-CO'^^'' reac- 
tion to be dominant), since within the mean- field picture 
reaction of the constant O*^' population with the linearly 
decreasing amount of CC^"*^ at the surface would give 
rise to a linear {\ — 9) dependence. With the surface 
kinetics of RuO2(110) dominated by ideal terrace sites 
and lateral interactions between O and/or CO at these 
terraces rather small Q , one would generally expect this 
mean-field analysis to be applicable, and - together with 
further experiments on which we comment below - the 
higher reactivity of the O™^ species was correspondingly 
questioned [31]. 

We scrutinize this conclusion with first-principles kMC 
simulations of the presented TPR data. This technique 
[2, M, combines the accurate description of the individual 
elementary processes and their interplay with an explicit 
account of the detailed spatial arrangement of the chem- 
icals at the surface. Here, we consider all non-correlated 
site and element specific desorption, diffusion and reac- 
tion events on the lattice spanned by br and cus sites, 
employing for each process first-principles rate constants 
computed by DFT and harmonic transition state the- 
ory [6|. This is thus exactly the same setup which was 
described in detail before and with which the experimen- 
tally measured steady-state activity of this surface could 
be successfully reproduced [61]. The TPR kMC simula- 
tions were carried out on (40 x 40) lattices (800 br, 800 
cus sites) with periodic boundary conditions, using the 
experimental heating ramp of 4.5 K/s. Modeling the dis- 
sociative adsorption of O2, surfaces with a defined ini- 
tial content of O*^"^ species were prepared by randomly 
filling pairs of neighboring free cus sites until the de- 
sired coverage was reached, and then filling all remaining 
empty cus sites with CO molecules, as well as filling also 
all bridge sites with O atoms. For such defined initial 
populations, the TPR spectra were simulated by mea- 
suring the amount of desorbing CO2 resolved in 10 K 
bins during the temperature ramp and averaging over 
several simulations with different random number seeds. 
Extensive tests verified that none of the specifics of the 
numerical setup, and in particular the finite lattice size, 
affected the quantities reported here. For the discussed 
total CO2 yields (integrated over the experimentally em- 
ployed temperature range of 170 - 600 K) sHghtly wrong 
desorption peak temperatures obtained as a consequence 
of the employed DFT energetics are also of lower con- 
cern, and we arrive at e.g. essentially identical results 
and conclusions when employing the refined DFT barri- 
ers reported in Ref. [7j. 
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FIG. 3: (Color online) Snapshots of the surface population in 
the first-principles TPR kMC simulation for the surface cov- 
ered initially with IML O''", 0.5 ML 0"="= and 0.5 ML CO'^"'. 
Shown is a schematic top view of one quarter of the simula- 
tion area of (40 x 40) surface sites, where the substrate bridge 
sites are marked by gray stripes and the cus sites by white 
stripes. Oxygen adatoms are drawn as lighter gray (red) cir- 
cles and adsorbed CO molecules as darker gray (blue) circles. 
Left panel: Initial population at T = 170K, Right panel: 
Population at T = 400 K. At this temperature, according to 
mean-field > 98 % of all initially present CO'^"'' species should 
already have been reacted off by the low-barrier 0'^"''-|-C0'^"'' 
reaction. Instead we find 40% of them to be still present, 
namely essentially all those that were not initially adsorbed 
immediately adjacent to a 0"^"° species. 



Performing our TPR kMC simulations for surfaces pre- 
pared with varying cus-site populations, we obtain ex- 
actly the same variation of the total CO2 yield with O*^"^ 
coverage as in the experiments, cf. Fig. [51 The linear 
dependence interpreted as evidence for a comparable re- 
activity of the O'"' species thus results despite the fact 
that the 0'^"''+C0'^"'^ reaction with a calculated barrier of 
0.9 eV [6] has a temperature-dependent rate constant in 
our simulations that is always orders of magnitude higher 
than the one of the competing O'^'^ + CO'^"'' reaction with 
a barrier of 1.2 eV [6[. Disentangling how much of this 
CO2 yield is produced by which of these two reaction 
mechanisms in Fig. [51 we see that the linear coverage de- 
pendence results in fact as the sum of a parabolic 6{l — 6) 
contribution due to 0'^"''-|-C0'^"'^ and a not quite linear 
(1 - 9) contribution due to O*"'' + C0™^ Rather than 
the qualitative functional form it is thus primarily the rel- 
ative weight of the two mechanisms that largely differs 
from the one deduced from the mean-field analysis. 

This is most apparent for the surface initially prepared 
with 0.5 ML O™" and 0.5 ML CO™", where mean-field 
would predict that essentially all CC^"** would be reacted 
off by the low-barrier O'^^'^-I-CO'^"'' reaction before the 
Qbr^(^Qcus j-gaction can efficiently set in at higher tem- 
peratures. Instead we obtain an almost equal contribu- 
tion of both reaction mechanisms to the total CO2 yield, 
cf. Fig. [21 which shows that the perfect mixing assump- 
tion behind mean-field is not fulfilled at this surface. In- 
triguingly, the corresponding inhomogeneity in the sur- 
face adlayer is here not created by lateral interactions 



(which are absent in our first-principles model [6[), but 
follows from the spatial arrangement of the active surface 
sites and diffusion limitations of the reactants. At the 
conditions of the TPR experiments the rows of strongly- 
bound O*^' species confine the reactive O'^^'' species to 
one-dimensional cus site trenches, cL Fig. [1] With the 
concomitant high diffusion barriers inside these trenches 
[g|, only a fraction of the CO'^"'' molecules can access 
the adsorbed O'^^^ atoms and react readily. This is il- 
lustrated in Fig. [31 which shows snapshots of the surface 
population at the beginning of the TPR ramp and at a 
temperature of T = 400 K. This temperature lies above 
the desorption peak temperature for the 0'^"''-|-C0'^"'' re- 
action as predicted by mean-field for a second order pro- 
cess and our temperature-dependent first-principles rate 
constant. According to mean- field essentially all of the 
initially present CO'^"*' species (> 98%) should therefore 
have already been reacted off. In our kMC simulations we 
find instead still 40 % of them to be present at the surface, 
namely essentially all those that did not have an imme- 
diately neighboring 0'^"'^ atom in the initially prepared 
ensemble, cf. Fig. [31 Due to the spatial arrangement 
of the two prominent active sites, all adsorbed CO*^"" 
molecules are on the contrary always chaperoned by two 
adjacent O*"^ neighbors, enabling a significant fraction of 
0'^''-|-C0'^^*' reactions despite the much higher reaction 
barrier. In this situation it is thus more the neighbor 
coordination, not the reactivity that determines the rel- 
ative weight of the two reaction mechanisms in the final 
CO2 yield, and this cannot be captured within a mean- 
field approach. Obviously, this geometric effect is spe- 
cific to the TPR experiments, since under steady-state 
conditions it would be efficiently counteracted by the re- 
plenishment of the cus population due to on-going dis- 
sociative O2 and unimolecular CO adsorption. We note, 
however, that even then a mean- field analysis fails qual- 
itatively, predominantly due to its inability to properly 
describe the site correlations introduced by the necessity 
to dissociatively adsorb O2 into two neighboring vacant 
surface sites [8[. 

When employing isotope labeling, the contributions of 
isqcus g^^ legbr gpgcies to the total CO2 yield can also 
be separated experimentally. Let us assume now that 
the corresponding two simulated curves in Fig. [21 were 
the data measured in a "computer experiment" . A rea- 
soning within the prevalent mean-field picture for TPR 
data analysis would have left us no alternative, but to 
conclude from the roughly similar contributions of both 
reaction mechanisms particularly aX — 0.5 ML that the 
0™"-f CO™" and O^^'-f-CO™" reaction had roughly equal 
rate constants. Already at this model surface of mod- 
erate complexity and without any lateral interactions, 
the mere neglection of the way the two prominent active 
sites are arranged at the surface alone would have thus 
led us to deduce rate constants that are several orders 
of magnitude wrong. Or, in terms of reaction barriers 



we would not have been able to resolve the 0.3 eV dif- 
ference in the barriers of the two mechanisms that was 
really the input to our "computer experiment" . The er- 
rors introduced by the mean-field assumption to the TPR 
data analysis are thus highly comparable to those ob- 
tained when using mean-field based rate theory to extract 
kinetic parameters from steady-state catalytic activity 
measurements j8|] . Such erroneous kinetic parameters are 
a source of confusion as exemplified by the conjectured vi- 
olation of the Br0nsted-Evans-Polanyi type relationships 
for the present RuO2(110) surface |3|, or - maybe even 
worse regarding theoretical work - this kind of data would 
then possibly be used to benchmark the quality of first- 
principles electronic structure calculations. 

Isotope labeling experiments of the just described 
type, i.e. separating the contributions of i^qcus ^^^ 
i6Qbr species, have in fact been performed and have in- 
deed shown the surprisingly low relative weight of the 
0™"-FC0™" reaction in the total CO2 yield [sf. The rel- 
ative magnitude and qualitative shape of the measured 
curve are similar to our simulation data shown in Fig. [21 
but its maximum is intriguingly shifted from 9 = 0.5 ML 
to ~ 0.2 ML. We are able to explain this finding as a con- 
sequence of an insufficient sample preparation. The ex- 
perimental data was recorded sequentially starting with 
the 6 = 0.8 ML TPR run. During this run a small amount 
of 18QCUS ^^j^ diffuse to bridge vacancies created by the 
legbr ^ CO™^ reaction. For the consecutive TPR run at 
the next lower coverage, the surface was only exposed to 
20 Langmuir ^SQa at 750 K (3,|3|. This achieves the in- 
tended annealing of remaining bridge vacancies, but does 
not remove the already present i^O^'' species. The en- 
suing amount of preexisting i^O^'' species at each run at 
lower coverages of i^qcus falsifies the obtained isotope- 
labeled CO2 yield and shifts the maximum of the curve, 
which we are fully able to reproduce when exactly fol- 
lowing the described experimental procedure [9| in our 
simulations [10|. 

In conclusion, we have shown that the detailed spatial 
arrangement of the active surface sites has a large effect 
on a recently reported set of TPR experiments study- 
ing the catalytic CO oxidation at RuO2(110). Even in 
the absence of lateral interactions at this surface, this 
nanostructure alone leads to inhomogeneities that cannot 
be grasped by prevalent mean-field data analysis proce- 
dures, which thus lead to misleading conclusions on the 
reactivity of the different surface oxygen species. When 
instead analyzed within an appropriate first-principles 
statistical mechanics framework the atomic-scale infor- 
mation provided by the present TPR data for RuO2(110) 
complies fully with the established picture of the func- 
tion of this model catalyst surface. This highlights the 
necessity of more refined data analyses when attempting 
the atomic-scale characterization of ever more complex 
surfaces exhibiting an ever wider variety of active sites. 
Correlations between the locations of these different ac- 



tive sites, i.e. their explicit arrangement at the surface, 
will then play an increasing role, making site accessibility 
a more and more important factor for the surface kinet- 
ics. The RuO2(110) surface studied here exhibits only 
two prominent active sites, arranged in simple alternat- 
ing rows. Yet, the mere neglection of this still quite triv- 
ial nanostructure leads mean-field TPR data analysis to 
extract kinetic parameters that are in error by several 
orders of magnitude and that do not even reflect the rel- 
ative reactivity of the different surface species correctly. 
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Insufficient sample preparation is also the presumed rea- 
son behind the observed desorption of isotope scrambled 
^*02 after exposing the stoichiometric RuO2(110) surface 
to ^®02 [S]- Preceding this experiment the surface was 
only annealed to 600K [3]. In our simulations we see that 
this is not sufflcient to desorb all O'^"'' species, since due 
to the severe diffusion limitations within the cus trenches 
there is always a fraction of isolated and at 600 K still 
practically immobile Q'^^" atoms remaining. These atoms 
can then only desorb at a much higher temperature either 
via recombination with a neighboring O'"^ species or be- 
cause diffusion is then facilitated and enables associative 
desorption of two remaining O'^"'' species. Such remnant 
16QCUS gpggjgg that are inadvertently present after the 
600 K annealing and therefore before saturating the cus 
sites with isotope- labeled 1*0'^"° provide a natural ex- 
planation for the measured fraction of isotope scrambled 
■^^02 in the TPD runs. The originally proposed explana- 
tion for these findings 3^ was an 0'^"''-0'"^ exchange dif- 
fusion process, which considering the low defect density 
at this surface would have to take place quite efficiently 
even at ideal terrace sites. Within our DFT setup [g| we 
compute a barrier above 2 eV for such a process, and can 
therefore rule out this original proposition. 



